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At sufficiently high temperature and density, quantum chromodynamics (QCD) 
predicts phase transition from the hadronic phase to the quark-gluon plasma phase. 
Lattice QCD is the most useful tool to investigate this critical phenomenon, which 
status is briefly reviewed. The usual problem in the Lagrangian formulation at 
finite density is either an incorrect continuum limit or its complex action and a 
premature onset of the transition as the chemical potential is raised. We show how 
the difficulties are overcome in our Hamiltonian approach. 

1 Introduction 

1.1 Motivation 

According to the standard model of cosmology, 12-15 billion years ago, the 
early universe underwent a series of drastic changes. Several microseconds 
after the big bang it was in a hot and dense quark-gluon plasma (QGP) 
state, where quarks and gluons were deconfined. Several minutes later, it 
changed to a phase, with the quarks and gluons confined inside the hadrons. 
Today it remains in the low temperature and low density hadronic phase. The 
ultimate goal of machines such as the Relativistic Heavy Ion Collider (RHIC) 
at BNL and the Large Hadron Collider (LHC) at CERN is to create the high 
temperature QGP phase and replay the birth and evolution of the universe. 
In February 2000, physicists at. CERN declared that they had experimentally 
found a new state of matter H, which had never seen before. Determining 
whether it is a QGP state requires further investigation from the future RHIC 
or LHC experiments. QGP may also exist at low temperature, in the core of 
very dense stars such as neutron stars. 

Quantum chromodynamics (QCD) is the fundamental theory of quarks 
and gluons. Figure 1 shows a recently proposedn two-flavor QCD phase di- 
agram on the temperature T and chemical potential fi plane. It assembles 
the information accumulated from lattice QCD, random matrix model, and 
other approximation methods. There is confinement and spontaneous chiral- 
symmetry breaking in the hadronic phase, which is separated from other 
phases by a chiral phase transition. The transition changes from second to 
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first order at a tricritical point. At high densities and low temperatures, a su- 
perconducting phase "2SC" was recently founda in which up and down quarks 
with two out of three colors form Cooper-pairs characterized by a diquark con- 
densate < >■ The transition between this 2SC phase and the QGP phase 
is likely first order. The chiral phase transition on the horizontal axis between 
the hadronic and 2SC phases is first order. 

T 



QUARK-GLUON PLASMA 




M 



Figure 1. Phase diagram of QCD with 2 flavors. 



A precise understanding of the QCD phase structure will provide valu- 
able information in the experimental search for the QGP. Lattice gauge the- 
ory (LGT) proposed by Wilson in 1974, is a first principle nonperturbative 
technique for the investigation of phase transitions. There are no free pa- 
rameters in LGT when the continuum limit is taken, in contrast with other 
nonperturbative techniques. Although the standard lattice Lagrangian Monte 
Carlo method works very well for QCD at finite temperature, it unfortunately 
breaks down at non-zero chemical potential (d ue to the so-called complex ac- 



tion problem). This is briefly reviewed in Sect. L2. In Sect. [2| we present our 
recently developed Hamiltonian approach fj to lattice QCD at finite density, 
by introducing the chemical potential in a natural way, which avoids those 
usual problems. 
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1.2 Present status 

In the continuum, the grand canonical partition function of QCD at finite 
temperature T and chemical potential fi is given by 

Z = Tr e -/3(H-MN) 5 /3= (] £bT )-1 ) (1) 

where ks is the Boltzmann constant, H is the Hamiltonian, and N is particle 
number operator 

N = [ d 3 x ^(x)tp(x). (2) 



The energy density of the system with free quarks is given by □ 

V Z V 0/3 W 

Going over to T — ► and chiral limit m — > 0, one gets the energy density 
(where the contribution of fi = is subtracted) 

Unfortunately, in Lagrangian LGT, a naive discretization of the chemical po- 
tential term does not lead to the correct continuum relation Eq.fl). Let us 
take the naive fermions as an example. The action reads 

a 3 - 

Sf = a 4 ^ mij}(x)ij}(x) + — tp(x)jk^(x + k) + a i iJ,}^ y ip^(x)ilj(x), (5) 

x x.k x 

where 7_fc = — 7^. In the chiral limit m — > and the continuum limit a — > 0, 
the subtracted energy density from this action is e su j oc (/i/a) 2 , i.e. becoming 
quadratically divergent and therefore it is inconsistent with the continuum 
result of Eq.(Q). This problem is not due to the species doubling of naive 
fermions, because the case of Kogut-Susskind fermions or Wilson fermions is 
similar. _ 

Hasenfratz and Karsch □ proposed to introducing the chemical potential 
exponentially 

3 3 

Sf = a 4 ^2 mip(x)ip(x) + \^( x )^i^( x + h ~ i>(x + j)ljiJK x ) 

X X j — 1 

a 3 

+ — [e^fcfr^te + 4) - e'^^x + l)^{x)] . (6) 
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The chemical potential can be introduced analogously for KS as well as for 
Wilson fermions. Such treatment of the chemical potential is numerically fea- 
sible in the quenched approximation (where the fermionic determinant det A 
is constraint to be 1, and quark loops are suppressed). However, there is 
evidence LI that the quenched approximation produces an unphysical onset 
of the critical chemical potential at the value fic — M n (m ^ 0)/2, which 
goes to zero in the chiral limit [M w (m ^ 0) is the pion mass at finite bare 
quark mass m. A finite bare quark mass has to be introduced in most of the 
numerical simulations]. This is in conflict with other theoretical predictions 
fic ~ Mffl /3 [M^ is the nucleon mass at fi = 0]. 

For full QCD, the fermionic degrees of freedom have to be integrated out, 
and in the measure occurs det A. For finite chemical potential det A from 
Eq. (H) and the effective fermionic action In det A become complex, due to 
Ay,x 7^ 75 Aj. 75, which renders numerical simulations extremely difficult. 
Much effort has been made to solve this notorious complex action problem: 

(1) The Glasgow group has suggested to treat det A as observable!! This 
method requires a very large number of configurations, in particular for /i rs 
He- Even on a very small lattice V = 4 4 , the computational costs exceed the 
current computer capacity H. Even in this case, the unphysical onset of nc 
still exists. Therefore, it is unclear whether the onset is an intrinsic problem 
of the proposal Eq. (^|). n 

(2) In the imaginary chemical potential method u det A becomes real, which 
works well for numerical simulations at high temperature and low density. 
But it might not work at low temperature and high density. 

(3) It has been proposed to utilize a special symmetry This is the only 
essful method in Lagrangian formulation, and attracts much attention 
, bat it works only for the SU(2) gauge group. _ 

(4) In Eil, the probability distribution function method t3 is applied to 2-color 
QCD (free of complex action problem), and very interesting results for the 
diquark condensate < ijnp > are obtained. 

2 Hamiltonian Approach 

2. 1 Free fermions at zero chemical potential 

The lattice Hamiltonian describing noninteracting Wilson fermions in d + 1 
dimensions at fi = is 

H = m{p{x)tp(x)+^2 (j^( x + k)+r (i/)(x) - ip(x + fe))) . (7) 

an x,k 
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The up and down components of aL&re coupled via the 7^ matrices, and one 
can use a unitary transformation E-3 

H' = exp(-iS) H exp(iS) (8) 

to decouple them. The operator S can be computed explicitly0 in momentum 
space: 




The physical vacuum state of H reads 

|n> = exp(iS)|0>, (10) 

where |0 > is the bare vacuum state defined as £|0) = 77 1 0) = 0. The parameter 
6 p in Eq. ([)[) is so determined that the vacuum energy Eq = (il\H\d) = 
(0|-ff'|0) is minimized. After the unitary transformation, the fermionic field 

tp can be simply expressed by up and down 2-spinors £ and as ip = y ^ 

and H' is diagonal 

h' = = E ( m + f - E sin2 fe a / 2 ) 

For Wilson fermions, in the continuum limit a — ► 0, for any finite momentum 
p, we have E n = —2N c Nf ^2 p a/to 2 + p 2 , giving the correct dispersion rela- 
tion. Here N c and Nf , respectively, are the number of colors and number of 
flavors. 

2.2 Free fermions at nonzero chemical potential 

Now we naturally introduce the chemical potential 

H^=H- fxN, (12) 

where H is given by Eq.(|^) and TV is given by Eq.(||). Let us define the state 
\n p ,n p ) by 

€ P \0p,n p ) = 0, ^,\0p,n p ) = |l p ,n p ), £ P \l P ,n p ) = \0 p ,n p ), ^,\l p ,n p ) = 0, 

Vp\n P ,Q P ) =0, i]l\n p ,0 p ) = \n p , l p ), r) p \n p , l p ) = \n p , P ), r)l\n p , l p ) = 0. 

The numbers n p and n p take the values or 1 due to the Pauli principle. By 
definition, the up and down components of the fermion field are decoupled. 
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Obviously, this is not an eigenstate of due to the non-diagonal form of H . 
However, they are eigenstates of H' which are related to by a unitary 
transformation 

Hi = exp(-iS) Hf, exp(iS) = H ' - fiN. (13) 
For the vacuum eigenstate of we make an ansatz of the following form 

= exp(iS)^2f np ,n p \n p ,n p ). (14) 
p 

S is given by Eq.Q, and H' is given by Eq. (fl))- The vacuum energy 
is E n = 2N c N f J2 P C npt n p [{A' p -fi)n p + (A' p + flfhp - A' p - p] . Here we 
have introduced the notation C np ,n p — fn p ,n p ' which have not yet been 
specified. For this purpose we use the stability condition of the vacuum. 
Because fi > 0, the vacuum energy increases with n p . This means the 
vacuum is unstable unless n p — 0. This simplifies the vacuum energy to 
Eq = 2N c Nf J2 P [Ci P {A'p — fi) — A' p — ju] , where we use the abbreviation 
C Up — Cn p ,o and the normalization condition. The dependence of C\ on 
the value of \i can be seen by inspection of the derivative dEo,/dC\ p = 
2N c Nf (A' p — /ij . For /i > A' p , the right-hand side is negative. Maximiz- 
ing C\ means minimizing the vacuum energy. Therefore, C\ = 1. On the 
other hand, for /i < A', the right-hand side is positive and for any C± p the 
vacuum is unstable. Therefore, C\ p — 0. We can summarize these properties 
by writing C\ p = G (ji — A' ) . Thus the vacuum expectation of H becomes 
(n\H\Q.) = 2N c N f J2 P (C lp A' p - A' p ). The subtracted energy density reads 

_ (n\H\n) - (n\H\n)i =0 _ ^ 

sub ~ N c NfN s ~4tt2- [W) 

Here N s is the number of spatial lattice sites. Thus we have proven that our 
Hamiltonian approach to free quarks at finite chemical potential leads to the 
correct continuum result for the vacuum energy density, Eq. (|J). From this 
relation, we can easily see that the free quark number density is proportional 
to /Li 3 . For naive fermions, in the continuum limit a = 0, there will be an extra 
factor of 2 d . 

2.3 Strong coupling QCD at nonzero chemical potential 

As is well known, lattice QCD at fi = confines quarks and spontaneously 
breaks chiral symmetry. For a sufficiently large chemical potential, this picture 

"There is a typo in Eqs. (2.24) and (2.25) in I: Eq there should be read as (Q\H\Q). The 
results and conlusions are the same. 



luo 'qcd: submitted to World Scientific on February 1, 2008 6 



T A 


1 


7j 


74 


75 


^7475 




le jjlj2ljllj2 




La 


1 


25 kJ - 1 


-1 


-1 


1 


1 - 26 kd 


25 k j - 1 


1 - 25 k j 



Table 1. T matrices and coefficients. 



may change. Here we set out to investigate finite .density QCD in the strong 
coupling regime 1/g 2 « 1. Following referenced, H' in Eq. @ is now 
replaced by 

x2 



l-(20 o rd + 



(20 o ) d 



4o 



32aiV r 



X] L ^^/i ( a; ) r AV , /2 (x)"0/ 2 (x + fyTAipft (x + k). 



(16) 



The additional four-fermion term is induced by gauge interactions with 
fermions. There d = 3 denotes the spatial dimension, fx, f% are flavor indices 
(summation over repeated indices is understood), 0q = 1 / '(4ma + g 2 ''C W) , and 
Cjv = (^V 2 — l)/(2-^Vc)- The matrices and their coefficients La are given 
in Tab.[l]. 

The unitary transformed Hamiltonian Eq. ( |l6|) can alsoJpe re-expressed 
in terms of the following pseudoscalar and vector operators Eil 

-.ift (1 - 74) 75^, 



n 



2v-w 
1 

' 2T-v 



■0 t (l-7 4 )7iV', 



(17) 



as I 



H> = 4 0) + Gx J2 nt (P)H(P) + E V l (P) y i (?) 
p \ i 

+ G 2 J] (nt(p)Ht(-p) + fc.c.) J] cosft-a 
p j 

+ G2 ( V j(p) V j(-P) + h - c ) cospya - 2 cospja 

P,3 \ j' ) 



where 



4° } = N f N s 



i-(2e yd 



(29 Q )d 



(18) 
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M M g 2 C N d(20 o f f g*C N (20 ) 2 d / t _ 

NfN s —v' — NfN s \vX — vo 

1 4a 1 l6aN r V 2 



d = 2 



l-(20 o ) 2 d 



(20 o ) d 



g 2 C N d(26 a ) 
AaN c 



G 2 



<7 2 CW (20 o )' 



8aiV c 



27V, 
"AT 



- J]] C np ,n p {Tip — Tip + 1) , 



2iVn 
(2JV C ) 



-J-jr- Gn p ,n p (Tip + Tip — 1) , 



t'2 



(2iV c ) 2 
A 7 . 



5]]^,% {n p + n p - \y 



Eq. ( |l8| ) can be diagonalized by a Bogoliubov transformation Q 
P{p) = coshu(p)a(p) + sinhu(p)a' (— p), 
Vj(p) = cosh. Vj(p)b(p) + sinh Vj (p)b^ (— p), 

-2G 2 



tanh 2u{p) 



tanh 2vj (p) 



E 



cospj-a, 



-2G 2 



cospjia — 2 cospj-a 



The Bogoliubov transformed Hamiltonian eventually becomes 

N 2 



II _ 77.(0) 



i?" = E, 



A r 



r Gi E 



1 - tanh 2u(p) - 1 



£g *E 



1 - tanh 2 2vj (p) - 1 



Gl E V 1 ~ tanh2 2u (pW (P)a(p) 



G i E V 1 - tanh2 2 «i(p)6}(p)6i(P)- 



(19) 
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Using the notation, normalization condition and arguments for the coefficient 



as in Sect. 2.2, we obtain n„ = and C\ 



i (0) 

l dyn 



Here to 



(0) 
dyn 



d/(ag 2 C]\[) is the dynamical quark mass at [i = 0. Therefore, in the chiral 
limit to = 0, the vacuum energy is 

Eq. 



2N c N f N s 



i (0) 

l dyn 



i (0) 

l dyn 



l (0) 
l dyn 



2^ Gl ^- 

p 

P-J 



1 - tanh 2 2u(p) - 1 



1 - tanh 2 2v 3 (p) - 1 



(20) 



According to the Feynman-Hellmann theorem, the chiral condensate is 
related to the ground state energy by 

1 ,. 8E n (to ^ 0) 



lim 



(0) 



NfN s m-To 9m 
where (ipip)( ' is the chiral condensate at /x = 



i-e /i 



i (0) 



(0) 



-27V, 



c 



1 



4fZ 



1 



N 



h 



Nf 



(21) 



(22) 



9 N ' 

and for d = 3, I x = 0.078354 ±2 x 10~ 6 , h = 0.235075 ±4 x 10~ 6 . According 
to Eq. ©, for /i < rnfy n , (ipip) — (iptp)^ ^ 0, i.e., chiral symmetry is 

spontaneously broken. For fi > m^ n , ('ipip) = 0, i.e., chiral symmetry is 
restored. Therefore, there is a first order chiral phase transition and the 
critical value of \i is given by 



(o) 

1 

dyn 



-2?^' ( 23 ) 

The critical chemical potential fic is equal to the dynamical quark mass, at 
/z = 0, which agrees with the result from an entirely different method E3 in 
Hamiltonian formulation. (The authors argued this was a second order phase 
transition, in contrast we clearly observe a first order transition). Our result 
is consistent with other theoretical predictions fic ~ because (see 

below) at fi = holds M^ 0) w im [ ^ n . 

We can compute now the quark number density in the chiral limit m = 0, 
which yields 

-1 dE n , (niE^CrcMaOlfi) 



2N c N f N s dfi 



1 



2N c N f N s 



l = QQi-Hc), (24) 
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This is consistent with the (3 = simulation results described in 0, and how- 
ever, is different from the large /i behavior in the continuum (i.e. the Stefann- 
Boltzmann law n q oc /i 3 for free quarks). It remains to be seen whether higher 
order 1/g 2 calculations will improve this behavior. 

The quark number susceptibility, standing for the response of the quark 
number density to infinitesimal changes in ji, is 

(in 

X 9 = ^ = <5(M-Mc). (25) 

Finally, let us look at some implications on the thermal mass spectrum 
of the pseudoscalar meson, vector meson and nucleon. The thermal mass is 
defined by M£ — (h\H — fiN\h) — Eq. For the pseudoscalar meson, in the 
chiral limit m = 0, 

for (i < fiQ, 



M* = G l Jl-t a nh'2u{p = 0)=i (0) r ^ (26) 

V \ Am dyn for M>MC 

Therefore, in the broken phase, the pseudoscalar is a Goldstone boson (M* oc 
yfm — > 0), and in the symmetric phase, it is no longer a Goldstone boson. For 
the vector meson, 

M v = G iy /l- tanh 2 2vj(p = 0) = ( ^fj ** (27) 

Am dyn for M > MC- 



where My ' 1 = A\/d — l/(a<7 2 Cjv) is the vector mass at fx = 0. Therefore, 
dM/dfJ, oc S(fi — He) for the pseudoscalar and vector mesons. To see the 
critical behavior at zero temperature, one should be very close to [ic- This 
behavior is consistent with that of the quark number density. To see whether 
the meson thermal masses depend on /z, higher order 1/g 2 corrections must 
be included. 

For the nucleon, we obtain the expected behavior 

M* N = M { ° } - 3^ (28) 
for (j, < fie, where w 3m^. This leads to Mj^ — at fx — He- 

3 Outlook 

In this paper, we have developed a Hamiltonian approach to lattice QCD at 
finite density. The chemical potential is introduced in a very natural way as in 
the continuum. It avoids the persisting problem in the Lagrangian approach 
of an incorrect continuum limit, or complex action or a premature onset of 
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the transition to nonzero quark density as [i is raised. The main result in 
the free case is given by Eq. (15), and those in the strong coupling regime 
are given by Eqs. (|2^)-(|2|). We have seen that the approach works well 
in the free case and also in the strong coupling regime. We predict that at 
strong coupling, the chiral transition is of first order, and the critical chemical 

potential jic ~ Afjy /3. 

We have not yet specified the nature of the chiral-symmetric phase for 
fx > lie- Is h a QGP phase or a color-superconduction phase 13 (see Fig. 1)? 
Up to now, there has been no first principle investigation of such a phase. 
Better understanding of it will affect our knowledge of the formation of the 
neutron star. 

At this point, we should mention that the results of Hamiltonian formula- 
tion and Lagrangian formulation La might be quite different. It is interesting 
to see whether they describe the same physics in the continuum. For pure 
gauge theory, within a Hamiltonian approach, we can extend to the intermedi- 
ate coupling and obtain meaningful results for the glueballs La. For fermions, 
the calculation is far from trivial. Recently we proposed a Monte Carlo tech- 
nique in the Hamiltonian formulation E9 for the purpose to do nonperturbative 
numerical simulations, by combining the virtues of the Monte Carlo algorithm 
with importance sampling and the Hamiltonian approach. We hope to apply 
it to QCD and with the aim to obtain useful information for RHIC and LHC 
physics and neutron star phenomenology. 



Acknowledgments 

We are grateful to the attendees of the workshop for useful discussions. 
X.Q.L. is supported by the National Science Fund for Distinguished Young 
Scholars (19825117), National Natural Science Foundation (10010310687), the 
Ministry of Education, the Doctoral Program of Higher Education and Guang- 
dong Provincial Natural Science Foundation (990212) of China. X.Q.L. and 
E.B.G. are supported by the Guangdong National Communication Ltd. H.K. 
has been supported by NSERC Canada. 



References 



1. http://press.web.cern.ch/. 

2. K. Rajagopal, Nucl. Phys. A661, 150 (1999), and references there in. 

3. E.B. Gregory, S. Guo, H. Kroger, X.Q. Luo, Phys. Rev. D62, 054508 
(2000). 

4. J.R. Li, Introduction to Quark Matter Theory, (Hunan Education Press, 
1989); J. Kapusta, Finite Temperature Field Theory, (Cambridge Uni- 



luo'qcd: submitted to World Scientific on February 1, 2008 11 



versity Press, Cambridge, England, 1989); B. Freedman and L. Mclerran, 
Phys. Rev. D16, 1130, (1977); 16 1147 (1977); 16 1169 (1977). 

5. P. Hasenfratz and F. Karsch, Phys. Lett. B125, 308 (1983). 

6. J. Kogut, M.P. Lombardo, and D. Sinclair, Phys. Rev. D51, 1282 (1995); 
I. Barbour, S. Morrison, E. Klepfish, J. Kogut and M. Lombardo, Phys. 
Rev. D56, 7063 (1997). 

7. I. Barbour, C. Davies, and Z. Sabeur, Phys. Lett. B215, 567 (1988). 

8. R. Aloisio, V. Azcoiti, G. DiCarlo, A. Galante, and A. Grillo, Phys. Lett. 
B453, 275 (1999). 

9. M. Alford, A. Kapustin and F. Wilczek, Phys. Rev. D59, 054502 (1999); 
M.P. Lombardo, Nucl. Phys. B(Proc. Suppl.)83. 375 (2000). 

10. S. Hands, J. Kogut, M.P. Lombardo and S. Morrison, Nucl. Phys. B 
558, 327 (1999). 

11. R. Aloisio, A. Galante, V. Azcoiti, G. DiCarlo, and A. Grillo, these pro- 



ceedings, |hep-lat/0007018| . 

12. Y. Liu, O. Miyamura, A. Nakamura and T. Takaishi, ticp-lat/0009009| . 

13. V. Azcoiti, V. Laliena, and X.Q. Luo, Phys. Lett. B354, 111 (1995). 

14. X.Q. Luo and Q. Chen, Phys. Rev. D46 (1992) 814; Q. Chen, and X. 
Luo, ibid. 42, 1293 (1990). 

15. J. Bjorken and S. Drell, Relativistic Quantum Mechanics (McGraw-Hill, 
New York); P. Roman, Advanced Quantum Theory (Addison- Wesley, New 
York). 

16. A. Le Yaouanc, L. Oliver, O. Pene, J. Raynal, M. Jarfi and O. Lazrak, 
Phys. Rev. D37, 3691 (1988). 

17. R. Aloisio, V. Azcoiti, G. DiCarlo, A. Galante and A. Grillo, Nucl. Phys. 
B564, 489 (2000). 

18. N. Bilic, K. Demeterfi, B. Petersson, Nucl. Phys. B377, 651 (1992). 

19. X.Q. Luo and Q. Chen, Mod. Phys. Lett. All, 2435 (1996); X.Q. Luo, 
Q. Chen, S. Guo, X. Fang and J. Liu, Nucl. Phys. B(Proc. Suppl.)53, 
243 (1997). 

20. H. Jirari, H. Kroger, X.Q. Luo and K. Moriarty, Phys. Lett. A258, 6 
(1999). X.Q. Luo, H. Jirari, H. Kroger, K. Moriarty, these proceedings. 



luo 'qcd: submitted to World Scientific on February 1, 2008 12 



